Quantum ratchet effects induced by terahertz radiation in GaN-based 

two-dimensional structures 
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I. INTRODUCTION 



In recent years, physics of quantum ratchets draw a 
growing attention. In mechanics, a ratchet is a device 
that is used to restrict motion in one direction while per- 
mitting it in another. In general, one means by a ratchet 
any sort of asymmetric potential, or a potential lacking a 
center of the spatial inversion. At the turn of 1980/90's 
it was understood that unbiased noncquilibrium noncen- 
trosymmetric systems can generate transport of particles, 
and this conception has roots in different fields of physics, 
chemistry and biology. In mechanical, electronic, opti- 
cal and biological systems, a particle, classical or quan- 
tum, charged or neutral, propagating in a periodic poten- 
tial with broken centrosymmctry and subjected to an ac 
force exhibits a net dc macroscopic flow. Examples are 
electrons in solids, Abrikosov vortices in superconductors 
and biological motor proteins^ One of implementations 
of the ratchet phenomenon is a Linear Photo-Galvanic 
Effect (LPGE). It represents a generation of a dc electric 
current under absorption of linearly polarized light in un- 
biased crystals of piezoelectric classes. Glass et al^ were 
the first to attribute a photoinduced current observed in 
ferroelectric LiNbOs to a novel photogalvanic effect and 
propose the first correct model for its microscopic qual- 
itative interpretation. A consistent quantitative theory 
of LPGE was developed by Belinicher and Sturman fol- 
lowed by other theorists, see Refs.yllj. 

Nonccntrosymmctric bulk semiconductors and het- 
erostructures are natural quantum ratchets and the fur- 
ther studies of LPGE in these systems allow one to eluci- 
date the whole problem of the ratchet effect. In particu- 
lar, a current of the LPGE consists of two contributions, 
ballistic and shift, which as a rule are comparable in the 
order of magnitude^ It is important to note that simi- 
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lar two contributions to the electric current, called the 
skew-scattering and side-jump currents, do exist in the 
anomalous Hall effect, see the recent review by Sinitsyn^ 
Another important point is that one should include into 
the light absorption process contributing to the ballis- 
tic current an additional simultaneous scattering from 
a static defect or a phonon, as compared to the light- 
induced optical transition leading to the shift current. 
In GaN-based structures first observation of LPGE has 
been reported quite recently il^ i 10 ! 11 ! 12 Here we present 
results of detailed experimental and theoretical investi- 
gation of LPGE in heterostructures based on GaN and 
its alloys with A1N. The commercial fabrication of blue 
and green LEDs have led to well established technological 
procedures of epitaxial preparation of these heterostruc- 
tures and initiated a great activity on investigations of 
their properties^ The photogalvanics serves as a solid 
bridge between transport and optics and, therefore, re- 
veals both transport and optical properties of the systems 
under study. 



II. SAMPLES AND EXPERIMENTAL 
METHODS 



The experiments arc carried out on (OOOl)-oriented 
wurtzite n,-GaN/Alo.3Gao.7N heterojunctions grown by 
MOCVD on C-plane sapphire substrates (for details of 
growth see Ref. 0). The thickness of the AlGaN layers 
was varied between 30 nm and 100 nm. An undoped 
33 nm thick GaN buffer layer grown under a pressure of 
40 Pa at temperature 550°C is followed by an undoped 
GaN layer (~ 2.5 /im) grown under 40 Pa at 1025°C; the 
undoped Alo.3Gao.7N barrier was grown under 6.7 Pa at 
1035°C. The mobility and density in the two-dimensional 
(2D) electron gas measured at room temperature are 
fi w 1200 cm 2 /Vs and N s « 10 13 cm -2 , respectively. 
To measure the photocurrent two pairs of contacts are 
centered at opposite sample edges with the connecting 
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FIG. 1: Photocurrent as a function of angle a measured at 
room temperature at oblique incidence (#o = 15°) in the 
a) longitudinal (J x ) and b) transverse (J y ) geometries. Pho- 
tocurrent is excited by linearly polarized radiation with wave- 
length A = 148 /jm and power P ~ 5 kW. Full lines are fits 
to Eqs. ||5J. To get agreement with experiments we used one 
fitting parameter Jo ex \ and introduced an offset J ff se t for 
the current J x detected in the direction of light propagation. 
The inset shows the experimental geometry. An additional 
inset in the lower panel displays the sample and the radiation 
electric field viewing from the source of radiation side. 



lines along the axes x || [1100] and y || [1120], see inset 
in Fig. [TJ The current generated by the light in the un- 
biased samples is measured via the voltage drop across a 
50 fl load resistor in a closed-circuit configuration. The 
voltage is recorded with a storage oscilloscope. 

The photocurrents were induced by indirect intrasub- 
band (Drudc-likc) optical transitions in the lowest size- 
quantized subband. The emission from a terahertz (THz) 
molecular laser optically pumped by a TEA CO2 laser- 
is used for the optical excitation. With NH3, D2O and 
CH3F as an active media we could obtain radiation pulses 
of duration ~100 ns with wavelengths A = 77, 90.5, 148, 
280, 385 and 496 fim. and a power P ~ 5 kW. Radia- 
tion in both the normal- and oblique-incidence geome- 
tries were applied with the angle of incidence 8q varying 
from —30° to +30°, 6q = corresponding to the normal 
incidence, and laser beam lying in the (xz) plane, see 
insets in Figs. [TJ H and g] 

Optically pumped molecular lasers emit linearly polar- 
ized radiation with the polarization plane determined by 
that of the pump light: the electric field vector of the 
THz radiation Ei generated by the molecular laser can 
be either parallel (like for the lines A = 90.5, 148 and 
280 fim) or orthogonal (lines 77, 385 and 496 /jm) to the 
polarization vector of the pump beam depending on the 
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FIG. 2: Photocurrent as a function of angle of incidence #0 
measured at room temperature in the transverse geometry 
and azimuth angle a = 45°. Photocurrent is excited by radi- 
ation with wavelength A = 148 fim and power P ~ 5 kW. Full 
line is the fit to Eqs. (|SJ with the same fitting parameter Jo 
as used in Fig. [TJ The inset shows the experimental geometry. 



angular momentum selection rules for the pump and the 
THz transitions in the active medial 

In the experiments the plane of polarization of the ra- 
diation incident on the sample was rotated applying A/2 
plates which enabled us to vary the azimuth angle a from 
0° to 180° covering all possible orientations of the electric 
field component in the (xy) plane. Hereafter the angle 
a = is chosen in such a way that the incident light po- 
larization is directed along the y || [1120] direction, sec 
inset in Fig.QJb). 

To investigate the photogalvanic effects we also use el- 
liptically polarized light. In this case the polarization 
of the laser beam is modified from linear to elliptical by 
means of crystal quartz A/4 plates. The counterclockwise 
rotation (viewing from the laser side) of the optical axis of 
the quarter- wave plate by the angle ip p results in the vari- 
ation of the radiation helicity as P C i rc = — sin2</? p . Here 
the angle tp p = is chosen for the position of the quarter- 
wave plate optical axis coinciding with the incoming laser 
polarization, in which case the linear polarization degree 
is given by fourth harmonics of ip p . 

A series of measurements is carried out making use 
of the frequency tunability and short pulse duration of 
the free electron laser "FELIX" at FOM-Rijnhuizcn in 
the Netherlands operated in the spectral range between 
70 /xm and 120 /im.— The output pulses of light from 
the FELIX were chosen to be m 6 ps long, separated by 
40 ns, in a train (or "macropulse" ) of 7 /xs duration. The 
macropulses had a repetition rate of 5 Hz. 



III. EXPERIMENTAL RESULTS 

Irradiating the (OOOl)-grown GaN/AlGaN heterojunc- 
tion by polarized light at oblique incidence, as sketched 
in the inset to Fig. QJa), causes a photocurrent signal 
measured across a contact pair. The width of the pho- 
tocurrent pulses is about 100 ns which corresponds to 
the THz laser pulse duration. The signal depends on 
the light polarization, and all characteristic polarization 
properties persist from 4.2 to 300 K. The experimental 
data presented below are obtained at room temperature. 
The effect is observed for all the wavelengths applied (be- 
tween 77 /xm and 496 /an). We use two geometries: the 
longitudinal geometry (J x in the Fig. [l}, in which the 
photocurrent is measured in the direction along in-plane 
component eii and the transverse geometry (J y in the 
Fig- U) j where the signal is detected in the direction nor- 
mal to the light propagation unit vector e [see the inset 
in Fig. [Ha)]. 

Figure [T] shows the dependence of the photocurrent 
on the azimuth angle a for both experimental geome- 
tries obtained at the positive incidence angle 9q = 15° 
and en || x. The polarization dependence of the cur- 
rent in the longitudinal geometry is well fitted by J x = 
Jq(1 — cos 2a) + Joffsct while for the transverse geometry 
we have J y = —Jo sin 2a. Note that the offset contribu- 
tion Joffset is observed in the longitudinal geometry only. 
Below, in Sec. IV. B, we will demonstrate that exactly 
these dependences follow from the theory. As a function 
of the incidence angle the photocurrent changes sign at 
9q w 0. This is demonstrated in Fig. [2] which shows the 
dependence J y (9o) for the fixed azimuth angle a = 45°. 
A similar dependence is also detected for the longitudi- 
nal current J. x . We note that, like the polarization de- 
pendent contribution to the photocurrent, the detected 
offset current in the longitudinal geometry J ffset inverts 
its direction when the incidence angle changes its sign. 

The wavelength dependence of the photocurrent ob- 
tained in transverse geometry at azimuth angle a = 135° 
and angle of incidence 0$ = 15° is shown in Fig. [3) The 
data are measured both on FELIX (dots) and on molec- 
ular laser (triangles). Fitting the data to the well known 
spectral behavior of the Drude absorptio n 16 i 17 
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we obtain that the spectral behavior of the photocur- 
rent can reasonably be described by this equation (see 
full line in Fig. [3|). For the fit we used r p as an ad- 
justable parameter and scaled the whole dependence by 
the ordinate. Analyzing this dependence we obtained 
that the momentum relaxation time controlling absorp- 
tion is about 0.05 ps. This value is twice shorter than 
the transport time (~ 0.1 ps) extracted from the mo- 
bility measurements at room temperature. We attribute 
the shortening of the momentum relaxation time by elec- 
tron gas heating due to absorption of the intense THz 
radiation. The reduction of the mobility due to heating 



GaN/AlGaN 
a =135° 




+ (05X^)2 



l" 




, >. 


I 15 




V^' c 


a 10 

O 




^ ' \ 



10 



r(K) 



100 



4 ^ 



o 2 



10-1 

50 100 200 300 400 

Wavelength, X (jim) 

FIG. 3: Spectral dependence of the transverse photocur- 
rent J y measured at room temperature at oblique incidence 
(#o = 15°) for azimuth angle a — 135°. The data are ob- 
tained with the free electron laser FELIX (dots) and with 
molecular optically pumped laser (triangles). Full line shows 
fit to Eq. Q. The fit is obtained using r v as an adjustable 
parameter and a scaling the whole dependence by the ordi- 
nate. Dashed line shows J y oc cu~ 2 for comparison. The inset 
in the left upper corner shows the temporal structure of the 
current in response to the radiation of FELIX. The inset in 
the right corner shows the temperature dependences of the 
carrier density and the mobility. 



stems from enhancement of scattering by phonons under 
increase of temperature. Using short 6 ps pulses of FE- 
LIX we observed that the response time is determined 
by the time resolution of our set-up but it is at most 
100 ps or shorter (see the inset in the left upper corner of 
Fig. [3]) . This fast response is typical for photogalvanics 
where the signal decay time is expected to be of the or- 
der of the momentum relaxation time^^ being in our 
samples at room temperature of the order of 0.1 ps. 

In addition to the signal excited by the radiation at 
oblique incident we also detected a small photocurrent 
at normal incidence. The polarization behavior of this 
signal is shown in Fig. 0J In this set-up photocurrent 
dependences can be well fitted by J x = J\ sin 2a and 
J y = J\COs2a (see also Sec. IV. A). Rotating the sam- 
ple around z axis we proved that the current direction is 
solely determined by the orientation of radiation electric 
field relative to the in-plane crystallographic axes. We 
note that this contribution is about one order of magni- 
tude smaller than that at oblique incidence even at small 
angles 8q w 15°. 

Besides investigations of photocurrent in response to 
linearly polarized radiation we also performed measure- 
ments under excitation with elliptically polarized light. 
Such experiments are of particular interest because such 
radiation, in partially circularly polarized light, has pre- 
viously been used for investigation of the circular pho- 
togalvanic effect (CPGE^i^ which coexist with the 
LPGE yielding beatings in <y9 p -dependcncc. The CPGE 
is characterized by a photocurrent whose direction is 
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FIG. 4: Photocurrent as a function of angle a measured at 
normal incidence (#o = 0) in a) x and b) y crystallographic 
directions. Photocurrent is excited by linearly polarized ra- 
diation with wavelength A = 148 fim and power P ~ 5 kW. 
Full lines are fits to Eqs. ((5}. To get agreement with exper- 
iments we used one fitting parameter J\ (X x'. The inset 
shows the experimental geometry. An additional inset in the 
lower panel displays the sample and the radiation electric field 
viewing from the source of radiation side. 



changed upon reversal of the radiation helicity. It has 
been observed in GaN-based structures demonstrating a 
substantial structural inversion asymmetry (SIA) in this 
wide band gap material ,L2jI2iii We have shown that SIA 
in GaN-based hetcrostructures results in a spin splitting 
of subbands in fc-space which was later confirmed by 
magneto-transport measurements! 18 i 19 i 20 i 21 The Rashba 
spin-splitting due to SIA, which is not expected in wide 
band semiconductors, in GaN/AlGaN heterostructures 
is caused by a large piezoelectric effect^ which yields a 
strong electric field at the GaN/AlGaN interface and a 
strong polarization induced doping effect *£• 

Figures OJa) and OJb) demonstrate the dependences of 
the photocurrent on the angle (p p for positive angle 8q. In 
Fig. EJa) transverse and longitudinal photocurrents are 
measured for the experimental set-up with Ei perpendic- 
ular to the incidence plane (s-polarization at <p p = 0) . We 
find that the polarization dependences in this experiment 
are well fitted by J x = Jq{\ — cos4</j p )/2 + J ff S ct and for 
transverse geometry by J y = —Jq sin4<^ p /2 + J 2 sin2(pp, 
see Sec. IV. B. 2 for theoretical justification. We note that 
we used for fitting the same values of Jo an d J ff S ct as in 
experiments with linearly polarized radiation. 

As we discussed above our laser can generate linearly 
polarized radiation either oriented along or perpendic- 
ularly to the polarization of the pump radiation. By 
that we change the position of the electric field of trans- 
formed beam relative to the optical axis of the quarter- 



wave plate, consequently changing the orientation of the 
ellipse as well as radiation helicity of the beam at the 
sample. In order to check how this transformation influ- 
ence results we provided an additional experiment using 
laser beam with again A = 148 /im and e | || x but for Ei 
parallel to the incidence plane (p-polarization at ip p = 0). 
Figures EJc),(d) show the data obtained in this geome- 
try. In contrast to the results presented in Figs.[HKa),(b) 
the photocurrent now shows another type of polarization 
dependences: J x = Jo(3 + cos4ip p )/2 + J Q Bset and for 
transverse geometry by J y = Josin4<^ p /2 + J2sin2ip p . 
The fact that contribution given by J2 did not change 
is not surprising because this term is due to radiation 



helicity which does not change: P c j rc 



sm2ip p . The 



changes in contribution induced by the degree of linear 
polarization are caused by the fact that the orientation 
of corresponding ellipses at the sample at a given ip p is 
different. 

The photocurrent polarization properties revealed ex- 
perimentally both under normal and oblique incidence 
are explained in the next Section where we propose a 
simple illustrative model to interpret the nature of the 
photocurrent under study and derive general phenomeno- 
logical equations. The measured wavelength and temper- 
ature dependences of the LPGE current are described in 
Sec. V where a microscopic theory is developed. 



IV. MODELS, PHENOMENOLOGY AND 

POLARIZATION DEPENDENCES OF THE 

PHOTOCURRENTS 

The appearance of the LPGE can be visualized by a 
simple one-dimensional model of randomly distributed 
but identically oriented wedges acting as asymmetric 
scattering centers! 3 ' 4 ' 14 This model showing the gen- 
eration of a photogalvanic current is relevant for 2D 
GaN/AlGaN structures excited by the radiation at nor- 
mal incidence. In Fig. [5] a wedge with a base oriented 
along x-direction is depicted which obviously does not 
possess a center of inversion. In equilibrium the veloci- 
ties of electrons are isotropically distributed. Application 
of an external in-plane alternating field E(t) = Esinuit 
adds an oscillatory motion along the electric field to the 
random thermal motion of the electrons. If the field 
points along the height of the wedges (E || y, or a = 0), 
then the scattering results in a carrier flow in the (—y) 
direction yielding an electric current j y > shown by 
up-arrow in Fig. [6^a). By that electron fluxes along x 
direction compensate each other and j x is absent. A vari- 
ation of the relative direction between the electric field 
and the orientation of the wedges changes the direction 
of the carrier flow resulting in a characteristic polariza- 
tion dependence. It may, e.g., reverse its direction, as it 
is shown in Fig. [BJb) for the field oriented along the base 
of the wedges, or rotate by 90° in the case of E directed 
at 45° to y axis. 

In order to describe the observed polarization and an- 
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FIG. 5: Photocurrent as a function of angle tp p measured at oblique incidence (#o = 15°) in the longitudinal [J x in a) and c)] 
and transverse [J y in b) and d)] geometries at wavelength A = 148 fim and power P ~ 5 kW. The ellipticity of the radiation is 
varied by passing linearly polarized laser radiation (Ei) through a quarter- wave plate (see insets in the left and right panels). 
Left panels show transverse and longitudinal photocurrents measured in experimental set-up with Ei perpendicular to the 
incidence plane (s-polarization at ip v = 0, see inset in the left panel). Photocurrents obtained for Ei parallel to the incidence 
plane (p-polarization at ip p = 0, see inset in the right panel) are shown in the right panel. Full lines are fits of the photocurrent 
to Eqs. (0 and (|10|) obtained correspondingly in the experimental set-up sketched in the inset in the left and right panels. The 
fits are obtained using the same values of Jo and Joffsct as in the experiments with linearly polarized radiation. The inset in 
the left panels shows the experimental geometry. On top the polarization ellipses corresponding to various phase angles ip p are 
plotted viewing from the source of radiation. 



gle of incidence dependences, we derive here phenomeno- 
logical equations for the photocurrents in 2D GaN-based 
structures and give a model for LPGE. The linear pho- 
togalvanic current density j is phenomcnologically de- 
scribed by the following expressions^ 



■LPGE 



= 2_, Xx^iE^E* + E V E*^ 

I"' 



(2) 



Here E is the electric field amplitude of the light wave, 
and x the third-rank tensor symmetric in the last two 
indices. The index A enumerates two in-plane coordi- 
nates x and y, while /i and v run over all three Cartesian 
coordinates. 

The LPGE is allowed only in piezoelectric crystal 
classes of noncentrosymmetric media where nonzero com- 
ponents of the tensor x do exist. The point symmetry 
of wurtzite 2D electron systems confined in the [0001] 
direction is C^ v which differs from the Cq v symmetry of 
the bulk GaN due to absence of translation along the 
growth axis z, i.e., because of presence of interfaces. The 
coordinate frame x,y,z chosen in the experimental set- 
up means that the plane (yz) || (1100) coincides with the 



mirror reflection planes a v contained in the C3„ point 
group. In this point group the tensor x nas two linearly- 
independent components for X =/= z: 



X = Xxx 



Xyyz 
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Xyx 



~Xyyy ■ V>) 



This means that the phenomenological equation ([2]) re- 
duces to 



LPGE 



X{E X E* Z } + X '{E X E*}, 



Jx 

3 L / GE - x{E y E* z } + X '{\E^-\E^) 



(4) 



where {E^E* V } = E^E* + E V E* 



A. Normal incidence 

In our experiments, a linear photogalvanic current is 
observed at normal incidence, see Fig. |4j It follows from 
Eqs. ([3]) that while at normal incidence the photocur- 
rent proportional to the constant x vanishes, because 
E z = 0, the contribution determined by the constant 
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FIG. 6: Model of current generation due to asymmetry of 
specular elastic scattering by wedges. The field E results in 
directed motion of carriers shown by dashed arrows. Due to 
asymmetric scattering a directed carrier flow and, therefore, 
electric current j are generated, a) and b) sketch two relative 
orientations of wedges and the electric field E together with 
resulting dc current j . 



X 1 is nonzero. The observation of this photocurrent is an 
important result demonstrating substantial difference in 
symmetry between wurtzite 2D GaN(0001)-based struc- 
tures with Cs v symmetry and zinc-blende GaAs- and 
InAs-bascd hetcrostructures having Ci v or Z?2d symme- 
try, where the LPGE is forbidden at normal light inci- 
dence. 

The polarization dependences for normal incidence fol- 
low from Eqs. (|4]) and are given by 



jx = x'(toE ) 2 sm2a, 
iy = x'(t £o) 2 cos2a, 



(5) 



where to = 2/(n ul + l) is the amplitude transmission coef- 
ficient for normal incidence, n^ is the refractive index of 
the medium (for GaN n w = 2.3). Both polarization de- 
pendences well describe our experimental data with one 
fitting parameter x' (see solid lines in Fig. 0|. This agree- 
ment clearly demonstrates the generation of LPGE cur- 
rent at normal incidence. We emphasize that the current 
direction depends on the orientation of the polarization 
plane in respect to the crystallographic axes x and y. 
In the case that contacts are arbitrary oriented relative 
to the crystallographic directions Eqs. ([5]) hold but the 
phase shift appears. This is considered in Appendix A. 



B. Oblique incidence 

1. Linear polarization 

Linear photogalvanic effect is also observed at oblique 
incidence. Moreover, comparison of Figs. [1] and 0] shows 
that even at a small incidence angle Qq the photocurrent 
is substantially (by the order of magnitude) larger than 
that at normal incidence. While at normal incidence the 
photogalvanic current is solely described by the constant 
x' at oblique incidence another term in Eqs. ([¥]) shows 



up. This term is determined by the second linearly inde- 
pendent constant X- The polarization dependence of this 
contribution in the geometry relevant to the experiment 
depicted in Fig.[TJ where the incidence plane is chosen to 
be the plane (xz) and the angle a is counted from y, is 
given by 



jx (a) = xEq t 2 p cos 6 sin 
Jv(<x) = -xEotptsS'm 



1 — cos 2a) , 
sin 2a . 



(6) 



Here 6 is the refraction angle related to the incidence an- 
gle 9q by sin (9 = sin#oAtaj and t s and t p are the Fresnel 
amplitude transmission coefficients from vacuum to the 
structure for the s- and p-polarized light, respectively^ 
These functions are shown in Fig. Q] by solid lines. To get 
agreement with experiments we used the above equations 
with one fitting parameter x and introduced an offset for 
the current J x detected in the direction of the light prop- 
agation. Registration of these characteristic polarization 
dependences proves observation of LPGE at oblique in- 
cidence. 

Comparing the results presented in Figs. [1] and U we 
can estimate of the ratio of different components of the 
tensor x f° r the studied structures. At 6 « 0.1 rad, 
t s f=a t p w ioj therefore from the ratio of the amplitudes in 
Figs.[T]andHequalto Ix'/x^l we get \ X '/x\ ~ 1CT 2 - This 
hierarchy of components of the third rank tensor XXftv 
in the systems of C% v symmetry is expectable^ however 
we demonstrate that the x'-related effect is observable. 
While the effect described by the constant x also exists in 
the bulk GaN as well as in the Coo V group the photocur- 
rent proportional to x' arises only due to the reduced 
symmetry C^ v of the system, i.e. due to size quantiza- 
tion. Therefore the latter contribution should increase at 
narrowing of 2D layer. 

Equations © show that the photogalvanic current 
should follow the dependence t s t p sin9, in particular, re- 
verse its direction upon inversion the angle of incidence. 
This behavior is indeed observed. Figure [5] shows LPGE 
current dependence on the incidence angle, 9q obtained 
for J y at a fixed linear polarization of radiation (a = 45°). 
We note that the offset photocurrent observed in the lon- 
gitudinal direction only is always directed antiparallcl to 
en. We ascribe the offset to the photon drag current 
which is linearly coupled to the photon momentum. This 
effect is out of scope of the present paper. 

In our experiments we also probe LPGE for the light 
directed along other crystallographic directions. Inves- 
tigating transverse and longitudinal currents being per- 
pendicular and parallel to the incidence plane, respec- 
tively, we observed that all polarization dependences are 
described by the same constant x independently of the 
in-plane direction of the light propagation. The reason 
of this fact is that the current contribution proportional 
to x is due to cone asymmetry being axial isotropic in 
the plane of our structures. Thus this photocurrent is to- 
tally determined by the light polarization state and prop- 
agation direction, and it is independent of the incidence 



plane orientation relative to the in-plane crystallographic 
axes. For linearly polarized light the photocurrent can be 
presented in the invariant form as j = 2xE\\E z where E« 
is the projection of the light polarization vector E onto 
the structure plane. It is seen that the LPGE current is 
always directed along E\\ . By introducing the azimuth 
angle (3 so that (3 = corresponds to the polarization 
perpendicular to the incidence plane (s-polarization) we 
obtain for transverse (j±) and longitudinal (Ju) currents 



j± = xEotpt s sin 2/3 (ey x z) , 
j'll = x^otp cos 9 (1 - cos 2/3)e|| 
where z is the unit vector along the normal. 

2. Elliptical polarization 



(7) 



and for the latter geometry (p-polarization at ip p = 0) 

i»(Vp) = 7;E$t 2 p cos6sin0 (3 + cos 4ip p ) , (10) 

jy(<Pp) = EfitptsSmO (- sin 4tp p + j sm2tp p ) . 

Here we again took into account that in our samples 
Ix'l ^ Ixl- Fits of experiment with the constant x i n ~ 
dcpcndcntly obtained from experiments with half-wave 
plates and 7 as a single adjustable parameter being the 
same for all curves yields a good agreement with the ex- 
periment. Here again for the longitudinal photocurrent 
the offset due to photon drag effect shows up and has 
the same value as that in the experiments with half-wave 
plates. 



Illumination of the structure by elliptically polarized 
radiation also results in generation of the LPGE current 
which is determined by the linear polarization degree of 
the light. However, elliptically polarized radiation is also 
characterized by the nonzero helicity. This results in the 
additional effect solely determined by the degree of cir- 
cular polarization, the circular photogalvanic effect.— ^^ 
This effect is phenomenologically described by 



■CPGE 

Jx 
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(8) 



where e is the light propagation unit vector and 7 the 
second-rank pseudotensor. In the systems of C^ v symme- 
try the tensor 7 has one linearly-independent component, 
namely ^ xy = —^f yx = 7- Thus, the CPGE current flows 
always perpendicularly to the incidence plane. 

In our experimental set-up a quarter-wave plate is ro- 
tated by the angle <p p between the laser light polarization 
vector, Ei, and the principal axis of the polarizer (see in- 
set to Fig. [5J. The total photogalvanic current is given by 
both second and fourth harmonics of this angle. While 
radiation handedness is solely determined by the angle 
ipp and does not depend on the relative orientation be- 
tween Ei and the plane of incidence the orientation of 
the ellipse is substantially different. In particular, for 
Ei perpendicular to the incidence plane at ip p — (s- 
polarization) and for Ei parallel to the incidence plane 
at ipp = (p-polarization) the orientation of correspond- 
ing ellipses on the sample differs by 90°. Therefore in 
the former case light is never p-polarized and in the lat- 
ter case s-polarization can not be achieved, which is quite 
different from the geometry of half-wave plate where all 
states of linear polarization can be obtained. Thus we de- 
rive two sets of equations describing experimental set-up 
ofFigs.EKa),(b) andl^c),(d). 

In the first geometry (s-polarization at tp p = 0) from 
Eqs. gl) and (g]) we get 



jx{<P P ) = -Elt 2 p cos 9 sin 6 (1 - cos 4^ 



(9) 
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V. MICROSCOPIC THEORY 

In order to describe the measured dependences of the 
photocurrent on wavelength and temperature we have 
developed a microscopic theory of LPGE. First we con- 
sider the contribution proportional to the coefficient x in 
Eqs. (4). This contribution does exist not only in systems 
of the C^ v or Cq v symmetry but is also allowed by the 
uniaxial point group Coo,,. Therefore, we can calculate 
the x contribution to j x p y GE in the uniaxial approxima- 
tion. Then we turn to calculation of the contribution 
proportional to x' i n Eqs. (4) and surviving under nor- 
mal incidence due to the reduced C3,, symmetry of the 
system. We note, that so far the LPGE in 2D systems has 
been treated theoretically for direct intcrsubband transi- 
tions only i 25 i 26 The mechanism of the LPGE under Drudc 
absorption differs strongly because scattering is unavoid- 
ably present in the intrasubband absorption process. 



A. Microscopic model for oblique incidence 

Microscopically the LPGE current consists of the shift 
and ballistic contributions which in general can be com- 
parable in order of magnitude. We start from the analy- 
sis of the shift contribution and then turn to the ballistic 
contribution. 



1. Shift contribution 

The shift contribution arises due to coordinate shifts 
of the free-carrier wave-packet center-of-mass by micro- 
scopic lengths in quantum transitions. The correspond- 
ing current is a product of the elementary charge e, the 
transition probability rate Wk'k and the shift vector Rk'k 
or, explicitly, 



sin Aip p + 7 sin 2(p p 



Jlpge = 2e Yl Wk ' k Rk ' 



(11) 



fc,fc' 



where the factor "two" accounts for spin degeneracy. The 
elementary shift is related to the transition matrix ele- 
ment Mk'k by^ 



Rk'k 



\M k , k \ 2 
= -(Vk + Vfc/)arg(M fc / fc ), 



(12) 



where arg(u) is the phase of the complex number u. The 
above two equations are general and valid for the shift 
photocurrent calculation in any frequency range^ii Here, 
for the first time, we apply them to consider the LPGE 
under the Drude absorption. Note that Eq. (12) was 
applied in the theory of anomalous Hall cffcct£>2£ where 
the coordinate shift at scattering is called a side-jump. 

The main contribution to the matrix element for indi- 
rect intrasubband optical transitions, 



M lfc',lfc = U lk>,lk 



Vw 



Ek' — Ek 



Vfe 



(13) 



comes from two-quantum processes ai and &2 with inter- 
mediate virtual states in the same el subband, see Fig. [7] 
Here the 2D electron kinetic energy Ek = Ti k 2 /2m, m 
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k' 



The processes ai and &2 in isolation from other pro- 
cesses make no contribution to JlpgE' because Afj fc/ lfe is 
independent of e z while the ^-related photocurrent in the 
phenomenological equations is proportional to e z , or 
because the microscopic expression in the square brack- 
ets in Eq. (|12p is real and hence the elementary shift is 
absent. The shift contribution becomes nonzero if two 
other indirect processes, bi and b2 in Fig. arc taken 
into account. The virtual states in these processes lie in 
the second size-quantized subband e2. The correspond- 
ing matrix element reads as 



eA 
M ik',ik = -ie z -^2 2 iA 2 i 
he 
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12 
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21 



V21 
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V21 



Ek' — E k 
(15) 
Here U21 = U12 is the intersubband scattering matrix el- 
ement, A21 is the energy spacing between the e2 and el 
subbands, and 021 is the intersubband matrix element of 
the coordinate z. The superscript _L indicates that the 
transitions (|15[) arc allowed in the polarization perpen- 
dicular to the interface plane. 

Substituting M k >k = ^\>k^ M k'k into Eq. (T2]), taking 



then into account the energy conservation equation (fTJ 
and assuming hw <C A21 we obtain for the elementary 
shift 



|M fc /, 



= -221 



Rk'k 

efiA 
me 



(16) 
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Note that the shift is different from zero only due to the 
wave vector dependence of denominators E k i — Ek and 
A21 + Ek> — Ek in the matrix elements (|13|) and (fT5"|) . 

The result for the shift photocurrent can be presented 
in the form 



FIG. 7: The quantum transitions responsible for the shift 
contribution to the photocurrent. Note that Drude absorption 
under oblique light incidence is caused by ai and a2 processes 
only. 

is the electron in-plane effective mass, Vk is the matrix 
element of electron-light interaction 

V k = e- k , 

mc 

u), A and e are the frequency, amplitude and unit po- 
larization vector of the light wave, Uik>,ik is the matrix 
element of intrasubband elastic scattering el, k — ► el, fc', 
the superscript |) indicates that the process is allowed for 
the polarization e containing an in- plane component e||. 
In the following, for simplicity we ignore the wave vector 
dependence of U\k\ik and replace it by a constant U\\. 
The probability rate Wk'k is given by Fermi's golden rule 
and expressed via the squared modulus of the transition 
matrix element and <5-function describing the energy con- 
servation 



Ek' — Ek 



huj 



(14) 
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Here the factor £ is given by-2£ 
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r (U11U21) 

where the angular brackets mean averaging over the 
distribution of static scatterers along z, the functional 
F[<fi(E)] turns any function <f>(E) into a real number as 
follows 



F[<j>(E)} 
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j[f{E)-f{E + hu>)]4>{E)dE I J f{E)dE , 



f{E) is the Fermi-Dirac distribution function, r/|| is the 
absorbance under normal incidence (6 = 0) given by2£ 



2nN, 
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F [2E + huj] 



(18) 



the momentum relaxation time t p = ft r / '(mNdiUfi)), N s 
and Nd are the 2D concentrations of electrons and scat- 
terers. 

Note that allowance for nonzero value of £ and, hence, 
of the LPGE current, presupposes the uniaxial asymme- 
try of the structure.— This corresponds to phenomenol- 
ogy of the LPGE which requires absence of an inversion 
center in the system. A value of £ is nonzero due to an 
asymmetrical shape of the electron envelope functions in 
the el and e2 subbands caused by the build-in electric 
field. An inhomogcncous asymmetrical distribution of 
scatterers along the z axis can also contribute to £. 



2. Ballistic contribution 

In order to obtain a nonzero ballistic photocurrent it 
is necessary to go beyond the approximation of single- 
impurity scattering and take into account simultaneous 
scattering by two impurities in one process. Such optical- 
absorption processes are schematically shown in Fig. [5] 
Let us introduce the antisymmetrical part, W k , k , of the 
fc — > fc' transition rate W k / k which changes sign under 
simultaneous inversion of the wave vectors fe and k' . It 
can be expressed in terms of Fermi's golden rule 



WS, 
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M k 'k\ a lf(E k ) - f(E k ,)] 5(E k , - E k 



fiuj) , 

(19) 

where |A/fc'fc| a is the antisymmetric part of the squared 
modulus of the optical two-impurity-mediated matrix ele- 
ment Mk'k- Then the ballistic contribution to the LPGE 
current is given by 
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T p (V k > - V k ) , 



(20) 



where v k = hk/m is the electron velocity. 

Here it is worth to discuss the dependence of the ballis- 
tic and shift contributions on the concentration of elastic 
scatterers Nd- For two-impurity scattering the proba- 
bility rate W k , k ex N%. In contrast, the single-impurity 
scattering rate W k ' k in Eq. (11) linearly depends on Nd- 
However, the ballistic current is controlled by the mo- 
mentum relaxation time, see the factor t p in the right- 
hand side of Eq. (20). Since t p oc N^ 1 , the ballistic and 
shift currents are both proportional to Nd- In this sense, 
there is an analogy between a pair of the two contribu- 
tions to the photocurrent and a pair of the side-jump and 
skew-scattering currents in the anomalous Hall effect £ 

In order to get the photocurrent proportional to both 
en and e z , we extract from |Mfc/fc| the interference term 
proportional to Re {M^M^*}, where M^ and M 1 - arc 
the contributions to the matrix element of the indirect 
optical transition proportional to e x ^ y and e z , respec- 
tively. To the lowest order in ?jw/A 2 i <C 1, one should in- 
clude into consideration six two-impurity scattering pro- 
cesses a 3 . . .ag contributing to M" and two processes b3, 
b 4 for M x depicted in Fig. [5] 



For example, below we present the contributions to M" 
from the processes as and a@, 
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where the prime means the scattering (shown by dotted 
lines in Fig. [5]) by a defect different from that which is 
related to the matrix element Un in these equations and 
to U21 in Eq. (p~5|) (shown by dashed lines). The Kro- 
necker symbols are obtained after averaging the matrix 
element over the in-plane distribution of the second scat- 
tcrcr. Note that averaging over its distribution in the z 
direction results in replacement of the square (U^) 2 by 
(f/ii). The energy denominators are rewritten by using 
the identity 

—^77 = PV- + inS{x) , 
x — iO x 

and, among three such terms in each contribution AfJ. 
(j = 3... 8), we retain the S- function in one of them and 
take the Cauchy principal values for two others. This 
additional (5-function, together with the energy conser- 
vation law (fLf| and the relation k + k^ = k' + fci , fixes 
absolute values of the wave vectors k, fc', fci and k-i- The 
integration over the azimuthal angle of the vector k\ ei- 
ther &2 is performed by using the principal value integral 
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where 8(x) is the Heaviside step- function. The calcu- 
lation shows that, among eight processes presented in 
Fig. [8j only two, namely, the processes as and ae, lead to 
nonzero ballistic current. The final result for this current 
reads 
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where 
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FIG. 8: Optical two-impurity-mediated processes included into consideration in order to obtain nonvanishing antisymmetric 
part of the scattering rate W^, k in Eq. (|19p . For the processes a3...6, the wave vectors satisfy the relation fc + ki = fc' + fci. 



The function L(x) goes to zero as x/4 at x — » 0, and 
tends to 1 at x — > oo. A sum of the shift and ballistic 
contributions given by Eqs. (17) and (20) represents the 
photocurrent phenomenologically described by the coef- 
ficient x m Eqs. (4). 



B. Microscopic model for normal incidence 

An ideal GaN/AIN interface has the point-group sym- 
metry C3,, which is a subgroup of the C|„ space group 
of a bulk wurtzite crystal. The group Cg„ is nonsym- 
morphic and contains, in addition to C3 and cr v , the el- 
ements {C6|r}, {(Td|x}, where r is the fractional lattice 
translation (0,0, c/2) and c is the lattice constant along 
the principal axis z. In the 2D structures, the screw-axis 
and glide-plane operations shift the interfaces along z not 
bringing the system into coincidence with itself and, as a 
result, the point-group symmetry reduces to C3,,. 

For the C61; symmetry of bulk wurtzite structure, char- 
acteristic of Ill-nitrides, the nonvanishing third order 
tensor components x\^v symmetrical in /1 and v are 

zzzi azii — Azyt/i aiiz — aizi — Ayyz — Xyzy — A- 
Therefore, forbidden are such bulk effects as the second- 
harmonic generation for excitation along z, a piezoelec- 
tric field induced in the (x, y) plane under the deforma- 
tion in this plane and the photocurrent in the (x, y) plane 
for the light propagating along z. The reduction from Cq v 
to C3,, removes these restrictions and four new nonzero 
components appear, \xx V = Xxyx = Xyxx = ~Xyyy = x'> 
see Eqs. (|4]). Note that it is the interface that is exclu- 
sively responsible for the appearance of new components. 

Here we propose a microscopic model of the LPGE 
under normal incidence. In this geometry virtual tran- 
sitions via other size-quantized electronic subbands are 
forbidden, therefore we should consider the light absorp- 
tion with intermediate states either in the same el elec- 
tronic subband or in the valence band. In the first case, 
see Fig. [7J the reduced C% v symmetry manifests itself 
in the scattering matrix element Uik',ik, and, instead of 
Eq. (HHJ), one should replace the elementary shift in the 



product |Mfe'fe| Rk'k by the scattering induced shift 
ju) _ M^ifc'.ifc^fe + Vfc/)C7i fa / )lfc ] 



R 



k'k 



\Uik',ik\ 



Note that the main contribution to U\k' ,ife dependent 
only on the difference (fc' — fc) does not result in the shift 
because the operator (V& + V&') nullifies such a term. 
Similar contribution to the LPGE current has been cal- 
culated for Drude absorption in the heavy-hole subband 
of a bulk p-GaAs^ where the shift appeared due to inter- 
ference of polar and deformation scattering mechanisms 
on optical phonons. 

For optical transitions going via T& valence-band 
states, the elementary shift can be related to asymmetry 
of the interband optical matrix elements in which case 
the asymmetry of the interband scattering matrix ele- 
ment is neglected. We consider this mechanism in more 
detail. Taking into account the symmetry considerations 
we can present the interband matrix elements of the op- 
erator e 11 • p (p is the momentum operator) between the 
conduction-band state Fi and the valence-band states Tq 
as 

e\\ ■ PcTi;v^ ex (fc) = '^Pe x + Q(e x k y + e y k x ) , (23) 
e\\ ■ p c r i;u ,r 6y (fc) = iVe v + Q{e x k x - e y k y ) , 

where V and Q are real constants. Multi-band deriva- 
tion of these contributions for GaN-based 2D struc- 
tures is given in Appendix [Bj The matrix element Q 
is nonzero due to the presence of the interface in the 
heterostructure. This approach is similar to descrip- 
tion of interface inversion asymmetry induced spin split- 
ting in the envelope-function method where ^-functional 
terms are introduced into the interband Hamiltonian, see 
Refs. I3lli32l . As a result, we expect a strong dependence 
of Q on the confinement size because Q is proportional to 
a product of the electron and hole envelope functions at 
the interface, see Eq. (|B3|) . We expand My p. and Rk'k 
in powers of small parameters huj/E g , E / ' E g <C 1. In the 
same nonvanishing order we should take into account k- 
dcpcndcncc of the matrix clement V '. This could be done 
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by replacing in Eq. (|23|) the constant V by the function 

Ek 
2E, 



V{k)=V Q [l--£- 



The interband scattering cT\k «-► vTe x k' and cTik <-> 
vTg y k' is described by the constant matrix elements U CVt . 
and U cv assumed to be independent of k and k' , with 

(ulj = (ul y ). 

In order to carry out a reasonable estimation of the 
coefficient \' m Eqs. ([4]) we have used a simple nonrcla- 
tivistic k p model coupling the cTi and vTq band states. 
In this model the valence band consists of two subbands, 
one with an infinite in-plane effective mass and the other 
with the hole in-plane effective mass coinciding with that 
for a conduction-band electron. Omitting the details of 
derivation we present the final result 
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(24) 



where F x = F[l] = f[f(E)-f(E+fiu>)]dE / J f(E)dE. 

o / o 

Experiment shows a noticeable magnitude of the LPGE 

current under normal incidence, see Fig. 2J 



VI. DISCUSSION 

In this Section we discuss the role of the contributions 
to the total current and analyze the frequency and tem- 
perature dependences of LPGE. An analytical estimation 
for the shift and ballistic photocurrents at Drude absorp- 
tion is readily available for the photon energies Tilo small 
as compared to the electron typical energy E which is 
relevant to our experiments performed at room temper- 
ature and applying THz radiation. In this case Eqs. \17\i 
and (|22[) reduce to 
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Since hw/E -C 1 the shift contribution dominates the 
LPGE. This estimation is confirmed by the study of the 
frequency dependence of LPGE. Indeed, while the shift 
photocurrent follows the frequency dependence of the ab- 
sorbancc which, for Drude processes, is given by Eq. |T]), 
the ballistic contribution has another frequency depen- 
dence because in this case Jlpge ^ r 7ii^-' 2 - The observed 
good agreement of the current and absorption frequency 
behavior (sec Fig. [3]) unambiguously proves the dominat- 
ing role of the shift contribution to the LPGE. We note, 
that such a spectral behavior of the LPGE has been pre- 
viously observed in bulk III-V semiconductors^ 

We note that even for the opposite limit of hco/E 3> 1 
the LPGE due to Drude absorption is still dominated 



by the shift contribution. In this case estimations yield 

.(shift) _ /o/q-v -(ball) 
JLPGE — \°/ Z JJLPGE- 

Now we analyze the temperature dependence. As fol- 
lows from Eq. (|17p . the temperature dependence of the 

LPGE current is given by Jlpge k £v\\(T). Taking into 
account temperature behavior of the Drude absorption 
[see Eq. ([15])] we get Jlpge <* £/ r p(^)- Here we again 
use the relation hu <C E ~ k^T. Variation of the tem- 
perature results in the change of the dominant scatter- 
ing mechanisms. Typically in the low-dimensional struc- 
tures by raising temperature we get a transition from 
the temperature-independent scattering rate by impuri- 
ties 1/Tp mp to the fast-increasing rate due to scattering 
by phonons 1/tS (T). As a result we have for the LPGE 
current 
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This expression demonstrates that despite the inverse 
mobility l/(i(T) ex l/rj mp + 1/tP h {T) is a monotonous 
function of temperature (see inset in Fig. [3]) and the 
asymmetry factors £ are temperature- independent, the 
LPGE current can have a minimum or even change the 
sign at temperature increase. This is due to the interplay 
between two contributions to the LPGE caused by dif- 
ferent scattering mechanisms characterized by substan- 
tially different values of £i mp and ^ p h possibly being even 
opposite in signs. Due to the fact that only one time 
is temperature dependent [t^ (T)] variation of temper- 
ature may result in the dominance of one or the other 
contribution to the LPGE current. 



VII. CONCLUSION 

In conclusion, we have studied the linear photogalvanic 
effect in wurtzite GaN(0001)-based heterojunctions un- 
der excitation by linearly and elliptically polarized light 
in the terahertz frequency range. Experimentally, the 
properties of LPGE have been revealed as a function 
of the incidence angle and the polarization state of ra- 
diation. A detailed point-group symmetry analysis of 
the polarization properties shows an agreement between 
the experiment data and the phenomcnological theory 
of LPGE. Phcnomenologically, the effect is described by 
two linearly independent coefficients x an d x', the former 
being allowed by the C3,,, Cq v and even Coo„ symmetries 
and the latter being nonzero only for the C3,, symmetry. 
An observation of the photocurrent under normal inci- 
dence determined by \' unambiguously demonstrates a 
reduction of the heterojunction symmetry from hexago- 
nal Cgi; to trigonal C3,,. However, beginning from very 
small oblique angles first contribution to the photocur- 
rent prevails over the second one. The frequency de- 
pendence of the measured LPGE photocurrent exhibits 
a Drude-like behavior. A microscopic theory of LPGE 
under the intrasubband (Drude) absorption in the lowest 
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el electron subband has been developed. While calcu- 
lating the ^-related photocurrent we have included vir- 
tual indirect optical transitions via the higher subband 
e2, estimated both the shift and ballistic contributions 
and showed that the first one exceeds the other by a 
large factor ex (E/Tilo) 2 . The dominating shift photocur- 
rent repeats a Lorentzian frequency dependence of the 
Drudc- absorption coefficient, in agreement with the ex- 
periment. The temperature dependence of the photocur- 
rent is governed by a product of the absorption coefficient 
and the dimensionless parameter of asymmetry. Micro- 
scopically, the contribution nonvanishing under normal 
incidence arises taking into account indirect optical tran- 
sitions via virtual electronic states in the valence band. 




FIG. 9: Schematic representation of the orientation of the 
crystallographic axes x, y relative to the directions x' , y' along 
which the current is measured. 



or, equivalcntly, 



Acknowledgments 

The financial support from the DFG and RFBR is 
gratefully acknowledged. E.L.I, thanks DFG for the 
Mcrkator professorship. Work of L. E.G. is also supported 
by "Dynasty" Foundation — ICFPM and President grant 
for young scientists. The high quality GaN samples 
were kindly provided by Hyun-Ick Cho and Jung-Hec 
Lee from Kyungpook National University, Korea. We 
gratefully acknowledge the support of the Stichting voor 
Fundamenteel Onderzoek der Materie (FOM) in provid- 
ing beam time on FELIX. 



jx> =x'(t E ) 2 sin (2a -3$), 
j y > = X '{t E ) 2 cos (2a -3*). 
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where a is counted from y': E(a = 0) 

Equations (|A1|) demonstrate that the azimuth angular 
dependence for arbitrary orientation of contacts in re- 
spect to the crystallographic directions remains the same 
but differs by the phase shift 3<&. 

Finally we note that the photocurrent determined by 
the constant x m Eqs. ^ is independent of the angle $ 
because it is totally defined by position of the incidence 
plane and polarization of light. Therefore Eqs. holds 
for any angle $. 



APPENDIX A: PHOTOCURRENTS IN THE 
ROTATED COORDINATE FRAME 



APPENDIX B: DERIVATION OF INTERBAND 
MATRIX ELEMENTS 



In the Section IV. A we get phenomcnological equations 
for LPGE at normal incidence. We obtained that the 
photocurrent direction is determined by the orientation 
of the electric vector with respect to the crystallographic 
axes x and y. Most convenient for this effect is to in- 
vestigate the current in the directions along and perpen- 
dicular to one of the mirror reflection planes of the C$ v 
point group. However, in contrast to III-V semiconductor 
heterostructures where well determined in-plane crystal- 
lographic orientation is naturally obtained by cleaving, 
in GaN/AlGaN heterojunctions grown on the sapphire 
substrate this may be not the case. Therefore we obtain 
here the LPGE for the arbitrary orientation of contacts 
in respect to the crystallographic directions. To describe 
this situation we introduce the angle <fr between the in- 
plane directions x' and y' along which the contacts are 
made and the crystallographic axes x and y, see Fig. |H1 
The symmetry considerations yield the expression for the 
combination of x' and y' components of the photocurrent 
which follows from Eqs. ^ for E z = 



Equation (f2"3"| can be obtained from the two-band elec- 
tron effective Hamiltonian 

Wcr 1) «,r 6 .(fc) = — (iVk x + Qk x k y ) , (Bl) 

m 



'HcT 1 ;v,T 6y {k) — — 

m 



Wk y + Q l -(kl-k 



by using the relation 



mo 



jx> + 



i.V = ix'C^*' - iE y ,)(E* x , - iB*,)exp(3i$) 



e-PoTi;ii,r 6 (fc) = -j~( e - V"fe)Hcr iru ,r 6 (fc) ■ 

We describe the symmetry reduction from C6u to C3,, 
by introducing a ^-functional perturbation V'5(z — Zjf), 
where za is the interface coordinate in the z axis and the 
operator V' transforms according to the representation 
Bi of the group Ce v with the basic function y 3 — 3x 2 y 
(here as before we assume that one of the three mirror- 
reflection planes in the C3,, group contains the axes y, z 
and is perpendicular to x). The perturbation V' mixes 
the Bloch state Tq with T$ and the Bloch state T^ with 
Ti. 



13 



In order to find the quadratic-in-fc correction in 
Eq. (|Bip we calculate the third-order correction to the 
Hamiltonian using the perturbation theory)^ 



<$>{ss') = 



L(£°-£°)(i?°-.E°) (£°-£°)(£°-£2) 



Here the indices c, v mean, respectively, the lowest con- 
duction and the highest valence bands of the symmetry 
Fi and Tq; s,s' are other conduction or valence bands 
different from c and v (the electron spectrum throughout 
the Brillouin zone for bulk GaN can be found, e.g., in 
Refs. [35ll36r ). 7i„ jra ' is the matrix of the first-order matrix 
elements of the k ■ p coupling or the interface mixing be- 
tween the states in the bands n and n'; E^ is the electron 
energy in the band n at the T point. The dimension of 
the matrix 7i n ,n' is N x N' where N, N' are degeneracies 
of the bands n and n' , respectively. Taking into account 
the bands sTq (s ^= v), sTj and s'T^ we can present 
the coefficient Q in Eq. (|B1|1 for a structure with a single 
interface in the following form 



V 



Q = VS(z - za) , 
h 



(B2) 



m n 



P'{R'V 2 - V X R) ${ss') - V 3 P"R ${s"s' 



The set of matrix elements is defined as follows: 



v s = (c,t 1 \v'\s",t 3 ), 



p' = -i^riifoia.rte) , p" = -i{s",r 3 \ Px \s',r 5x ) 



R= -i(s',r5a;|p x |u,r6a;) , R' = —i(s', T 5x \p x \s, T Qx ) . 

For the band sT 5 we choose the basis \s, T^x), 
\s,T 5y ) which transforms under operations of the C3,, 
group as the coordinates x and y, i.e., as func- 
tions \s',T 6x ), \s',T 6y ). While deriving Eq. (JB2J) we 
took into account the relation between matrix ele- 
ments imposed by the symmetry, e.g., (c,ri\p y \s,Te y ) = 
(c,Ti\p x \s,T 6x ), (c,T 5y \p x \s,T 6x ) = -(c,T 5y \p y \s,Te y ) 
= (c,T 5x \p x \s,T 6x ) etc. 

It follows then that, for a heterojunction, the inter- 
band optical matrix elements are indeed given by Eq. |23|) 
where 

V = P J f e (z)f h (z)dz , Q = Vf e (zu)Mzu) , (B3) 

P = -i(c,Ti\p x \v,T 6x ) = -i(c,ri|p w |u,r 6 i,), and f e (z), 
fh(z) are the electron and hole envelope functions. 



r, 



(s,r 6iC |t>'| s ',r 5:c ) , v 2 = (v,t 6x \v'\s',t 5x ) 
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